From 1437 runs
Filtering the parameters sets to match:
output years are 2005 to 2025 (climate repeats after 8/2018)
This is test for eucalyptus. Note that am not using the BAAD carbon ratio fractions here.
From sensitivity analysis (sobol.html), the top parameters that affect pre-drought summer LAI include:
## Warning: `funs()` is deprecated as of dplyr 0.8.0.
## Please use a list of either functions or lambdas:
##
## # Simple named list:
## list(mean = mean, median = median)
##
## # Auto named with `tibble::lst()`:
## tibble::lst(mean, median)
##
## # Using lambdas
## list(~ mean(., trim = .2), ~ median(., na.rm = TRUE))
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_warnings()` to see where this warning was generated.
We’ll first look at the spread of the model output.
The first pass of filtering results in 9291 parameter sets.
See how there’s two humps in the density plt? the idea is to get the second hump, that actually contains reasonable values of LAI. We can try that by filtering out the starting values in 2011 with plantc and LAI remote sensing estimates.
The second pass of filtering results in 4754 parameter sets.
let’s explore parameter values.
## $title
## [1] "leaf_turnover"
##
## attr(,"class")
## [1] "labels"
## `summarise()` ungrouping output (override with `.groups` argument)
Create new .def file for each species able to find
Test new .def file with spinup and run over drought years 2010-2018
Compare LAI / leafC with David’s index data over drought period 2011-2014-2017